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GENERAL THEORY OF AERODYNAMIC INSTABILITY AND THE MECHANISM OF 

FLUTTER 


tty Theodore Theodorsen 


SUMMARY 

The aerodynamic forces on an oscillating airfoil or 
airfoil-aileron combination of three independent degrees 
of f reedom have been determined . The problem resolves 
itself into the solution of certain definite integrals , which 
have been identified as Bessel functions of the first and 
second kind arid of zero and first order. The theory , 
being based on potential flow and the Kutta condition, 
is fundamentally equivalent to the conventional wing- 
section theory relating to the steady case. 

The air forces being known , the mechanism of aerody- 
namic instability has been analyzed in detail. An exact 
solution, involving potential flow and the adoption of the 
Kutta condition , has been arrived at. The solution is of 
a simple form and is expressed by means of an auxiliary 
parameter k . The mathematical treatment also provides 
a convenient cyclic arrangement permitting a uniform 
treatment of all subcases of two degrees of freedom. The 

flutter velocity , defined as the air velocity at which flutter 
starts, and which is treated as the unknown quantity, is 
determined as a function of a certain ratio of the fre- 
quencies in the separate degrees of freedom for any magni- 
tudes and combinations of the airfoil-aileron parameters . 

For those interested solely or particularly in the numeri- 
cal solutions Appendix I has been prepared. The rou- 
tine procedure in solving numerical examples is put 
down detached from the theoretical background of the 
paper. It first is necessary to determine a certain number 
of constants pertaining to the case , then t o* perform a few 
routine calculations as indicated. The result is readily 
obtained in the form of a plot of flutter velocity against 
frequency for any values of the other parameters chosen. 
The numerical work of calculating the constants is sim- 
plified by referring to a number of tables , which are in- 
cluded in Appendix I. A num ber of illustrative examples 
and experimental results are given in Appendix II. 

INTRODUCTION 

It lias been known that a wing or wing-aileron struc- 
turally restrained to a certain position of equilibrium 
becomes unstable under certain conditions. At least 
two degrees of freedom are required to create a con- 
dition of instability, as it can be shown that vibrations 


of a single degree of freedom would be damped out by 
the air forces. The air forces, defined as the forces due 
to the air pressure acting on the wing or wing-aileron 
in an arbitrary oscillatory motion of several degrees of 
freedom, are in this paper treated on the basis of the 
theory of nonstationary potential flow. A wing- 
section theory and, by analogy, a wing theory shall be 
thus developed that applies to the case of oscillatory 
motion, not only of the wing as a whole but .also to 
that of an aileron. It is of considerable importance 
that large oscillations may be neglected; in fact, only 
in fin i tel v small oscillations about the position of 
equilibrium need be considered. Large oscillations 
are of no interest since the sole attempt is to specify 
one or more conditions of instability. Indeed, no 
particular type or shape of airfoil shall be of concern, 
the treatment being restricted to primary effects. The 
differential equations for the several degrees of freedom 
will be put down. Each of these equations contains a 
statement regarding the equilibrium of a system of 
forces. The forces are of three kinds: (1) The inertia 
forces, (2) the restraining forces, and (3) the air forces. 

There is presumably no necessity of solving a general 
case of damped or divergent motion, but only the 
border case of a pure sinusoidal motion, applying to the 
case of unstable equilibrium. This restriction is par- 
ticularly important as the expressions for the air force 
developed for oscillatory motion can thus be employed. 
Imagine a case that is unstable to a very slight degree; 
the amplitudes will then increase very slowly and the 
expressions developed for the air forces will be appli- 
cable. It is of interest simply to know under what 
circumstances this condition may obtain and cases in 
which the amplitudes are decreasing or increasing at a 
finite rate need not be treated or specified. Although 
it is possible to treat the latter cases, they are of no 
concern in the present problem. Nor is the internal 
or solid friction of the structure of primary concern. 
The fortunate situation exists that the effect of the 
solid friction is favorable. Knowledge is desired con- 
cerning the condition as existing in the absence of the 
internal friction, as this case constitutes a sort of lower 
limit, which it is not always desirable to exceed. 
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Owing to the rather extensive field covered In the 
paper it has been considered necessary to omit many 
elementary proofs, it being left to the reader to verify 
certain specific statements. In the first part of the 
paper, the velocity potentials due to the flow around 
the airfoil-aileron are developed. These potentials 
are treated in two classes: The noneirculating flow 
potentials, and those due to the surface of discon- 
tinuity behind the wing, referred to as “circulatory” 
potentials. The magnitude of the circulation for an 
oscillating wing-aileron is determined next. The 

v 


VELOCITY POTENTIALS, FORCES, AND MOMENTS OF 
THE NONCIR CL'LATOR Y FLOW 

We shad proceed to calculate the various velocity 
potentials due to position and velocity of the individ- 
ual parts in the whole of the wing-aileron system. 
Let us temporarily represent tin* wing by a circle (fig. 
1). The potential of a source * at the origin is given 
by 

<p=r Io s y 2 ) 

4 7 r 

For a source e at (/ tt ?/i) on the circle 



/ 

•W s o 


Simjkk 1. Conformal represeniulimi of t lie wing profile by a circle. 

forces and moments acting on the airfoil are then 
obtained by integration. In the latter part of the 
paper the differential equal ions of motion are put 
down and the particular and important ease of un- 
stable equilibrium is treated in detail. The solution 
of the problem of determining the flutter speed is 
finally given in the form of an equation expressing a 
relationship between the various parameters. The 
three subcases of two degrees of freedom are treated 
in detail. 

The paper proposes to disclose the basic nature of 
the mechanism of flutter, leaving modifications of the 
primary results by secondary effects for future investi- 
gations. 1 Such secondary olTccIs are: The effects of a 
finite span, of section shape, of deviations from poten- 
tial flow, including also modifications of results to 
include twisting and bending of actual wing sections 
instead of pure torsion and deflection as considered in 
this paper. 

The supplementary experimental work included in 
Appendix II similarly refers to well-defined elementary 
cases, the wing employed being of a large aspect ratio, 
nondeformable, and given definite degrees of freedom 
by a supporting mechanism, with external springs 
maintaining the equilibrium positions of wing or wing- 
aileron. The experimental work was carried on 
largely to verify the general shape of and the approxi- 
mate magnitudes involved in the theoretically pre- 
dicted response characteristics. As the present report 
is limited to the mathematical aspects of the flutter 
problem, specific recommendations in regard to prac- 
tical applications are not given in this paper. 

1 Tin* filled uf intprnul 'riclion is in sr>nip cusps essential; this subject will be 
contained m a subsequent paper. 


,p= tw !o " K x “- r i)"+ ?/i) 2 } 

Putting a double source 2t ut and a double 

negative source —2c at (/ lt - t/, ) we obtain for the flow 
around the circle 


v ~ 


2tt 


I r f £ 1 /llljj ?/- ?/l)~ 

(7 - A)' --(?/ + yt) 2 


The function <p on the circle gives directly the sur- 
face potential of a straight line pq, the projection of the 
circle on the horizontal diameter. (See fig. 1.) In 
this case ?/= v 1 — and ip is a function of x only. 

We shall need the integrals: 


JM 


(r-s,y+(r/-wy 


Or -si ) J + 0/ + 2/d 


^ d/i =2 (./• - c) log A — 2 A 1 — x 2 cos 


and 


r\ (r r, c * 

Jc (/ s x y-\- 


; iy ?/i ¥ , 


where 


, , (j'j — Cpd^i = — V l — C'v 1 J- 

oos^cOr — 2c) y I —x 2 - 1 {i - c) log N 

Y ^ 1 fj - v'l --T y! 


The location of the center of gravity of the wing- 
aileron x a is measured from <i , tin* coordinate of the 
axis of rotation (fig. 2); x$ the location of the center 



FuU'ke 2.— Parameters* ol the airfmJ-aileroo combioaiion. 


of gravity of the aileron is measured from c, the coordi- 
nate of the hinge; and r a and r A are the radii of gyration 
of the wing-aileron referred to a, and of the aileron 
referred, to the hinge. The quantities and r$ are 
“reduced” values, as defined later in the paper. The 
quantities a , x a , c, and s 0 are positive toward the rear 
(right ), h is the vertical coordinate of the axis of rota- 
tion at a with respect to a fixed reference frame and is 
positive downward. The angles a and j3 are positive 
clockwise (right-hand turn). The wind velocity r is to 
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the right and horizontal. The angle (of attack) a I 
refers to the direction of r, the aileron angle 0 refers to 
the undeflected position and not to the wind direction. 
The quantities r a and r & always occur as squares. 
Observe that the leading edge is located at —1, the 
trailing edge at +1. The quantities a> c, x a , r a , 
and r &i which are repeatedly used in the following 
treatment, are all dimensionless with the half chord b 
as reference unit. 

The efTect of a flap bent down at an angle j8 (see fig. 
2) is seen to give rise to a function <p obtained by sub- 
stituting — r/3 b for c; hence 

^ = ^fVl — x 2 cos-'c- (x- c) log N] 

To obtain the effect of the flap going down at an 
angular velocity /3, we put, r - — (x, — c)/36 2 and get 

(pp = [ V 1 ~ C 2 V 1 — X 2 + COS“ l c(x — 2 C) \\ ~ X~ 

- (x c) 2 log N\ 

To obtain the effect of an angle a of the entire nir- ( 
foil, we put c= — 1 in the expression for <p fiJ hence 
<p a = vab y 1 — x 2 

To depict the airfoil in downward motion with a veloc- 
ity h ( + down), we need only introduce - instead of a. 

- v 

Thus . _ 

n = M>V 1 -x l 

Finally, to describe a rotation around point a at an 
angular velocity a, we notice that this motion may be 
taken to consist of a rotation around the leading edge 
c = — 1 at an angular velocity a plus a vertical motion 
with a velocity — a(l -\-a)b. Then 

ifa = 7r(x 4- 2) -y/l — x J — a(\ 4- ( l)b 2 y 1 x 2 
Ztt 

<l^X 1 x 2 

The following tables give in succession the velocity , 
potentials and a set of integrals 2 with associated con- 
stants, which we will need in the calculation of the air 
forces and moments. 

VELOCITY POTENTIALS 

v? a = vab yl — x 2 
(ph = kb v /i — x' 2 

ipi = ab 2 (^x a^x \ x 1 

<P 0 * * v&b{ V l — x 2 cos" l c — (x — c) log N\ 

7T 

<p$ = }- /&> 2 [ V 1 ~ c 2 y 1 - x 2 4- (x ■ • 2c) V 1 “ x 2 cos" 'c 

-7 r 

-- (x ~c) 2 log A 7 ] 

where ^ = l C x - Vl — V I - r 

x — c 

2 Some of the more <itl11eult inlegnil evaluations are triven in Appemliv Ul. 


INTEGRALS 


| *p a i\x 

j y?;>dx = 2^ iTa 


I ^ ^ li <U=2 raT 

r +i h i 

J_j ^*dx 2 tn 


j: 


<^dx ab 2 7 9 


<p,dx - ab' 2 ™ 


+1 b 


n b 

J, ~ 2r V ^ Tr> 

- - 2 -&T 2 J_ 1 Vddx= - y/37’, 


£ <p tt (x-c) 


— c)dx = — ^vaT, 


= — ^VaCir 


£ iph(x - c)dx = - \bT\ 
j' iPa(x-c)dx = ab‘ 2 T l3 
£ <p 3 (x - c)dx = - 

/> 2 • 

j y>i(x c)dx = 2^/87 3 


r+i 

I V?a(^ _ c)dx 

f+1 6 . 

I I ^(j* -c)dx= - 2 Acjt 

r +l 

I c)dx = a 6 2 7 T M ir 

C+ 1 b 

I 2 V P J '» 


CONSTANTS 

Ti - — g Vl ~c-’(2 + c 2 ) 4- c cos _1 c 

7V c(l -c 2 ) y i - c 2 (l -Fc 2 )cos“ 1 c4-c(cos"’c) 2 

T ? ^~ - f (ros _1 e) 2 4- “C Vf ~ c * cos _I c(7 4- 2c 2 ) 

-g(l-c 2 ) (oc 2 4 4) 

r J\ - cos _1 c I c V 1 — c 2 

T-, — (1 — c 2 ) — (cos - 'c) 2 + 2c yl — c 2 cos“‘c 

7V -7; 

74 — ^ + c 2 ^ cos" 1 c + £ c V 1 "" c 2 (7 4- 2c 2 ) 


7 g x 1 c 2 (2c 2 4 1) It cos 1 c 

r, = \ [|( V 1 - + a VI = l ( ■ - P + a 7’,) 

where p— -4 ( V 1 — c’) 1 

T|0 “ Vl cos" 1 c 

T u = cos" 1 c ( 1 - 2c) + v i — c 2 (2 — c) 

T 12 = Vl “ C 2 (2 4- c) - cos" 1 c (2c 4 1 ) 

7’„-5[-7’,-(e-a)7’,l 

rr _ 1 , I , 

10^2 (lC 

* FORCES AND MOMENTS 


The velocity potentials being known, w T e are able to 
calculate local pressures and by integration to obtain 
the forces and moments acting on the airfoil and 
aileron. 
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Employing the extended Bernoulli Theorem for un- 
steady flow, the local pressure is, except for a constant 


Ph = 


'(? + to 


where w is the local velocity and <p the velocity poten- 
tial at the point. Substituting w = v + ^ we obtain 

ultimately for the pressure difference between the 
upper and lower surface at x 




where v is the constant velocity of the fluid relative to 
the airfoil at infinity. Putting down the integrals for 
the force on the entire airfoil, the moment on the flap 



-AT 


Figure 3.— Conformal representation of the wing profile with reference to the 
circulatory flow. 

around the hinge, and the moment on the entire air- 
foil, we obtain by means of partial integrations 


P - - 
A/a= - 
M a - - 


r + 1 

2 pb I v?dx 

- b 2 ^ ${x — c) d/ + 2 pvb | <pdx 

r\ i , r+i 

2 pb- <p(x ~ c) dr -t 2 pvb ydx 

r+i 

- 2pb 2 I <p(c — a) dr 

Or, on introducing the individual velocity potentials 
from page o, 

P = — pb 2 [t’7ra + tt/i — bird tit — v — b ]\$\ (1) 

M, = - p6 3 [ - rT,& - T,h + 2T a ba - \ vTj - 1 T 3 6/3 J 


+ pvb 2 




l\h + ‘2TM 


X vTS~-TM 

7T T 


*] 


= — pb 2 ^7 V ,2 a — (2 + Ti) bra + 27 1 13 6 2 d + ^ T b r 2 p 

+ T 2 -± T 2 ) bvfi - 1 b l T& 4 7>A - Ti&jt J (1 1 ) 

A/ a - pb 2 ^ — irV 2 a + + fl 2 ^ P'a + / l\ ft + { 1\ — 7\ 

- (c — a) 7\ } b&o -t- { — 7 1 ? - (c a ) " [ b 2 $ 

— baich — irvh J (III) 


VELOCITY POTENTIALS, FORCES, AND MOMENTS 
OF THE CIRCULATORY FLOW 

In the following we shall determine the velocity 
potentials and associated forces and moments due to a 
surface of discontinuity of strength U extending along 
the positive x axis from the wing to infinity. The 
velocity potential of the flow around the circle (fig. 3) 
resulting from the vortex element —AT at ( X 0t 0) is 

„ r = fi - tnn-'-jjj 


AF * - 

— ;v tan 


X- 

(rb x -) r 


X 2 


( Xo+ i) 


\X+ P+ 1 


where (A”, Y) are the coordinates of the variable 
and Xq is the coordinate of — AT on the r axis. 

Introducing X Q + -i- = 2r 0 

or X Q — r 0 + -\/r 0 2 — 1 on the r axis 
and X=x and Y = \ 1 — x 2 on the circle 
the equation becomes 

AF 4 yT^Vro 2 ” 17 ! 

2ir tan ' 1-x/o 

This expression^ gives the clockwise circulation 
around the airfoil due to the element — AF at r 0 « 


We have: p 


~ p \i>t + dx) 


But, since the element — AF will now be regarded as 
moving to the right relative to the airfoil with a 
velocity v 

c V 


dt dx 0 


Hence, p =■ 
Further 




AF da: 


\J"o 1 


X 

;^o)V] 

l + 0 -J 2 ) (y 0 ; 


_ , /oyl -x 2 
x 2 (l-xx 0 ) 2 


1) 


, v 


/V- 1 


(1 -XX 0 ) 2 

1 


and 


Vl -P (x 0 “ x) 

1 0*0 


2t d^_ — 3 ( \~xx 0 ) y ar 0 2 — 1 

A I’’ At > 1 J / 1 _ 


+ W-1 


(1 - XXq) 2 


AF dj 0 


Bv addition: 


1 + 




(l-^o) 2 


= A 


l-~r 2 1 


y'V-i (x 0 — x) 

<V d^^_AF XqYx 

too" 7l^V?7- 
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To obtain the force on the aileron, we need the 
integral 

Zo+3* 


jr(S5*&)- £f 


~=cos l x+ 
Va*o 2 - 1 V*i 


V*o 2 - lVl-ar 

VT^ 7 


dx 


jz£1 1 

o 1 


-oos- c + ^ /l - c H 
1 Vxo 2 -iJ 


- A ■*; u “ a A< 

Thus, for the force on the aileron 


a/ j c ,= -pt’^(- / 4==co S -'c+ 1 4=-r Vi-c 2 ) 

T-VVZt) 1 V^O 1 / 


A/ J cl = — pvb 


Integrated, with AT = L^dxo 


3ci= _^[ ( 


cos 1 c~* y 1 — c 


*>r 


X(, 


V-Co 2 - 

fVPIH 


£/dx 0 


for c — 1 we obtain the expression for P, the force 

on the whole airfoil 


P = — pr6 


■jr 


3" n 


V*a 2 ~ 1 


Z7dx 0 


(IV) 


Since £/ is considered stationary with respect to the 
fluid elements 

U=i(vt — Xo) 

where / is the time since the beginning of the motion. 
U is thus a function of the distance from the location 
of the first vortex element or, referred to a system 
moving with the fluid, U is stationary in value. 

Similarly we obtain for the moment on the aileron 

f/|, + |^x b= Ar n .(x^cUxo +£ , 

Jc \i>x faj 2*Jc 

AI’ i r r, 5 , 

= ~ ^ vi7-TU oVl_ ^ + -T- — cvl ~ x 

+ (^-Ioc)cOS _1 xJ 


+ 2~. c )Vl“C 2 


,au 1 r, 

,y j 2 7 1 V^O 

- x \Zn ~ 1 L 

+ ^1 “2^oC^COS _1 cJ 

AT /ft / ,7 -5 

= + ■ ,-^r vv\i-c 2 ~cc 

2jrL\/o— 1 




os ’c) 

1 

VI -V 2 ) J 


Finally 

- cos-'f (c + ^)J + \ -] (cos-'c - cVl -c J ) ] 

Putting AT= £/<Lr 0 and integrating 

*--■£[>-^0+1) 

-eos-> C (c+|)JJ“^^'dx 0 

+ (cos-'c -cV(l-c 2 )|J VvFT L ' d/0 ] (V) 


Further, for the moment on the entire airfoil around a 


C + Y*£4.*e\ 

J— i U* W 

+ G~*" rt ) 


, , , Ar 

(x-a)dx— — ~ 


2,r Vx,r 




\ , n +i *r i 

/l \ 

)cos'*j_ r + 2„ Vxo >_l 


V / 


1 

2 _ ’ ;r o a 

and A M a = — pvb 2 AT 2 ~ - 

y x 0 — l 

Integrated, this becomes 



THE MAGNITUDE OF THE CIRCULATION 

The magnitude of the circulation is determined by 
the Kutta condition, which requires that no infinite 
velocities exist at the trailing edge, 
or, at x— 1 

! v* + v>d+- + <p$) = finite 

Introducing the values of <p a , etc. from page 5 and 

, 

from page 6 gives the important relation: 
l =va + h + b( i - a)d 

-ttJi yx 0 - 1 / 

+ —>'0+b^0 (VII) 

-7T 2 7T 

This relation must be satisfied to comply w T ith the 
Kutta condition, which states that the flow shall leave 
the airfoil at the trailing edge. 

It is observed that the relation reduces to that of the 
Kutta condition for stationary flow on putting j 0 = co , 
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and in subsequence omitting the variable parameters 
d, $, and i. 

Let us write 


_l r i z±} 

- T .Jl \x«-l 


L djr t) I’a f h I b 




+ 7 ’ u 't-(3d - Q 

7 r It 


Introduced in (i V) 

r-~L'dr 0 

r= - 2irpvbQ - 


from (V) 

Mp = — l 2pvb 2 

r j> 

Ji \ JV - 1 


rv 


/ Jo + 1 

Jo -1 


LVl/n 


J 


/ jp^ 1 

A Jo - 1 


■i( cos ' 


Ui. lx, 


c c V 1 — c 


Udxa 


( i+ j)) 
s ) 


X 



and from (VI) 


d/ a — — 2 irpr 6 2 


Introducing 


we obtain finally 


f 


Jo 


i 

ft A,. 1 - 1 

r " 


Vjq-~ l 

/ Jo + 1 

A 1 ' Jo - 1 

l dJo 


r ^ dj 0 


Udx o 


/j()^ 1 ITI 

\r 7 “d <Ia 


P — — 2pr6ir(7<^ 


(VI II) 


A/ fl - -2pi’6 2 ^Vl -c^l + 0-cos' 1 c(c + t))c 
+ ‘ (cos- 1 c - c v'l -T !! )Jv : prb’(T„C- T,)Q (IX) 

•V, - 27rpr6-’j^'i + l )o~ Q (X) 


where $ is given above and C^C(k) will be treated in 
the following section. 


VALUE OF THE FUNCTION C (k) 


Put (/== cv’b ('- “Ld 

where s ~ vt (s^co), the distance from thej£r$< vortex 
element to the airfoil, and k a positive constant deter- 
mining the wave length, 
then 


( \k) 


j; 

i 


Jo 

■ ( 

i J(T“ i 


T Oiljo 


Jo + 1 -ffc* ] 
-■■-—.e ** 2 odj 0 
A Jo* — 1 


(XI) 


These integrals are known, see next part, formulas 
(XIV) — (XVII) and we obtain 3 


<■(*)=- 


— j - 

. ? / __ n \ r 


,7T 

iu J 1 


1 71 7 

' * .>«A> 


-./, + Q’ . 

— {J i + ) o) + i ()*,—«/()) 


( t/, • i I' i) [ _(«/]_+_ ) o) v_{ ) ,) ] 

( r A ) n ) * -r ( I i — t/y ) 2 

J i (</ 1 ( i o) 4 - y i ( y i ~ e/ p) 

(^i “i* y o) z ’b ( y i _ t/o ) 2 


. } i (</i 4- y o) </i ( y i Ju) 

i*j\ y o ) 2 _ b ( y i — </o)* 




where 


a t = ^ 1 ^ 1 ^ i o) 4 ) i(y i t/p) 
(«A 4- Jo ) 2 4 (),— t/p ) 2 

/■ _ y i y o 1 e/ n/p 

(t/i + ) o ) 2 4- ( ) ! - t/p ) 2 


(XII) 

(XIII) 


These functions, which are of fundamental import- 
ance in the theory of the oscillating airfoil are given 

graphically against the argument ^ in figure 4 . 


SOLUTION OF THE DEFINITE INTEL KALS IN V BY MEANS OF BESSEL 
FUNCTIONS 


We have 


I< n 



e 


— z cosh l 


cosh nt d< 


(Formula (34), p. 51— Ciray, Mathews 
Sc MaeRobert: Treatise on Bessel 
Functions. London, 1922) 

where 

i n* 

Kn (t) = € 2 G n (it) 


(Eq. (28), sec. 3, p. 23, same reference) 

and 

On (j) = - Y„ (j) + [log 2- 7 + (j) 

but 

U u) = | Y, U) + (log 2 - 7 ) J n (x) 


(where Y n (j) is from N. Nielsen: 
Ilandbuch der Theorie der Cylinder- 
funktionen. Leipzig, 1904). 

1 This may al*i> In? «x [tressed in Hhnkel functions. C’- 


296 



GENERAL THEORY OF AERODYNAMIC INSTABILITY AND THE MECHANISM OF FI. OTTER 


G„ (r) — ? [ Y„ (x) iJ„(r) ] 


f 'o> r ® 

. J, v?-i 


0 el 1 j*il j: 


x , .T , , fN f - cos Avdx , ■ f"sin kxdx 

~r,U(k) + *r, X (*) - J, -,p=rr + * J, ^ f 

1US, 

pco M = _r Yo{lc) (XIV) 

txv) 

,/l -yr 1 

F urther, 

/' w 1 J*llx 

K,( ik) - | #-’* cosl1 ' costikl/ j i 

i'6', (A-j (k) 

= I - (cos 4 / sin A*/)!!/ 

Ji aV 1 

ms, 

(XVI) 

Ji \x 2 -- 1 

J X S 1 1 1 k X{ I X IT * r , » , y t t [ | \ 

— t= 2 — , r; 5 i (k) (XV II) 

i V-f "' 1 

TOTAL AERODYNAMIC FORCES AND MOMENTS 


TOTAL FORCE 


From equations (J) and (VIII) we obtain 

P - pb 2 (vna 4- nil ~ nb(la ~ vT 4 & — 1\ bfi) 

— Li 7 t pvbCJ j va 4 /<- + 6^ “ <^a4- ^;7’ lo r0 

+ h L T n& } (X 


TOTAL MOMBNTS 


From equations (II) and (IX) we obtain similarly 
M„ - pb[ [ 2 Ty-T, I ^(tf- *-)}»6a I'lTuVa 

+ -v-per, T,T W ) — vbf}T,T n - \rjM 

TT - 7T TV 

- - prb 2 TJ!\ va -I- k + iQ - a)a 

+ — 7'iiii^ I b^ r T„(t\ (X I X 


From equations (I l L) and (X) 

Af a p 6 ^ 7 r^ ^ i + rt 2 ^ci 

4- (7 t 4 H- Tn,)v 2 fi 
4 ^7\ - 7 T 8 - (c a) 7 4 + 

-- ^7 7 4- (c — rt) 7\ - anbh .J 

4- 2p?7) 2 7r^« 4- va 4 h + b(\ - 

4- ~Ti 0 vf} + b~- T u ji \ (XX) 

TT ~ 7 T J 


DIFFERENTIAL EQUATIONS OF MOTION 

Expressing the equilibrium of the moments about a 
of the entire airfoil, of the moments on the aileron 
about c, and of the vertical forces, we obtain, respec- 
tively, the following three equations: 
a : - I a a - - b(c - a) S$ - SJt - aC a 4- Ah = 0 

—I — Ipa — b(c — a) aSp — hS$ — fiCp 4- Alp = 0 
h : - Ml- ixS a ~ PS , - hC h 4 P = 0 

Rearranged: 

a : ah + fi {h 4 Vbifi — a) Sp) 4- XaS** 4 aC y a ~ Ah “ 0 

/3 : a(/^4- 6(c — a)Sp) 4- 4 hS$ 4- fiCp — il/^ -- 0 

h 1 4* fiSp 4~ hAl + hO n — P ■— 0 

The constants are defined as follows: 

Pj mass of air per unit of volume. 

b, half chord of wing. 

AI, mass of wing per unit of length. 

S a ,Sfi static moments of wing (in slugs-feet) per 

unit length of wing-aileron and aileron, 
respectively. The former is referred to 
the axis a; the latter, to the hinge c. 
I a ,h, moments of inertia per unit length of 

wing-aileron and aileron about a and c, 
respectively. 

C at torsional stiffness of wing around a, cor- 

responding to unit length. 

Cp, torsional stiffness of aileron around c, cor- 

responding to unit length. 

C nt stiffness of wing in deflection, correspond- 

ing to unit length. 

DEFINITION OF PARAMETERS USED IN EQUATIONS 



the ratio of the mass of a cylinder of air of 
a diameter equal to the chord of the 
wing to the mass of the wing, both taken 
for equal length along span. 
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■V& 

S a 


tlu* radius of gyration divided by b. 


rj- 


the center of gravity distance of the wing j ^ 


from a , divided by b. 


A i -r- j the frequency of torsional vibration 


iC a 

w a ^ -%/ — ’ 

around a. 

= -y/ reduced radius of gyration of aileron 
divided by 6, that is, the radius at 


y^} reduced center of gravity distance from c. 

I C~ 

frequency of torsional vibration of aileron 

around c. 

1C 


r» = 


— -y/ ^ frequency of wing in deflection. 

FINAL EQUATIONS IN NONDIMENSION AL FORM 

On introducing the quantities Af a , Mp, and P y 
which the entire mass of the airfoil | replacing and In from page 5, and i educing to 
would have to be concentrated to give nondimensional form, we obtain the following system 
the moment of inertia of the aileron la. of equations: 

w & [ r « + *(| + ° ! ) +& l K (l -a ) +o< M" 2 ~^[ r ' ,2 ' l ‘ ( 2 ~ *) 7< ] 

■f^i(A+A„) + 

(Bj a^r' 


+ (c-a)xa-K 




(c- 


- a)^_ +&(p-2\- l r <)ll + ii(r 0 2 - Va) - . yh 0AAij* 


T C a . 1 r* ,, r rrrp r \ , }•/ 1 T \1 , A« I'ClA-jf ra , h , /I \ . T m v T tl 

+ ? p K { h ~ T * T '°>] + b {**--* r ) i [tt " & T L ~ ■ a ) a " t b &+ T ^\ - 0 


(C) a + a-j k + $ (xe — ~ A*^ — A* j + A (1 + *) % + hj 

vC(k)\"v<x h ( 1 ? n a! _ n 

6 [v + S + U " V" + T If * + 2^_T 0 


c.i 

A/& 


SOLUTION OF EQUATIONS 

As mentioned in the introduction, we shall only have 
to specify the conditions under which an unstable 
equilibrium may exist, no general solution being 
needed. We shall therefore introduce the variables at 
once as sine functions of the distance s or, in complex 

form with ^ as an auxiliary parameter, giving the 

ratio of the wave length to 2ir times the half chord b: 

ifc4- 

a — a 0 e 0 

i ( k r + 

0=00* K b J 

and h = k Q e ^ h 

where s is the distance from the airfoil to the first 

dLs 

vortex element, ^ = v t and and <p ^ are phase angles 
of 0 and h with respect to a. 


Having introduced these quantities in our system of 


equations, we shall divide through 



We observe that the velocity v is then contained in 
only one term of each equation. We shall consider 
this term containing v as the unknown parameter ilX. 
To distinguish terms containing X we shall employ a 
bar; terms without bars do not contain X. 

We shall resort to the following notation, taking care 
to retain a perfectly cyclic arrangement. Let the 
letter A refer to the coefficients in the first equation 
not containing C{k) or A", B to similar coefficients 
■ of the second equation, and C to those in the third 
equation. Let the first subscript a refer to the first 
| variable a, the subscript 0 to the second, and It to the 
third. Let the second subscripts 1, 2, 3 refer to the 
second derivative, the first derivative, and the argu- 
ment of each variable, respectively. *4 a i thus refers 
to the coefficient in the first equation associated with 
the second derivative of a and not containing C{k) or 
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X; C\z to the constant in the third equation attached to 
h, etc. These coefficients 4 are as follows: 


A °^ r ; + (l +a *) 

' i "- = (2 — a ) 

d o3 = 0 

An = — — — + (c — a) (— — — 1 

K 7T \ K T / 

A»-\(T<+T W ) 


Af,\ — 



Am — 0 
A h3 ~ 0 


2»..-^-£+ ( c-a)(3 

K 7T \ K 

r.) 



5 o3 ' 0 


(=^l) 


The solution of the instability problem as contained 
in the system of three equations A, B, and C is given 
by the vanishing of a third-order determinant of com- 
plex numbers representing the coefficients. The solu- 
tion of particular subcases of two degrees of freedom 
is given by the minors involving the particular co- 
efficients. We shall denote the case torsion-aileron 
(a, fi) as case 3, aileron-deflection (0, h) as case 2, and 
deflection-torsion (A, a) as case 1. The determinant 
form of the solution is given in the major case and in 
the three possible subcases, respectively, by: 


D= 


^ao r~j~ Ra0~\~ Ha0i Rah +»/ ah 

Rba ~f~* i-ffitr j R j %i bp , Rbh +U bh 

Rca~\~Hca) 5 c 0 + i/ c 0, Rch~t~ H ch 


= 0 


and 


M ch = 
Af aa 


Raa + i ‘^aoi Ra$-\-HaP J 

Rba-^ril ben Rb&^ril b&\ 
j Rb&-\-iIbflf Rbh^ilbh 
\Rc$ + iIc,, RchA-ilch 
Rchflil Chi Rea +u 

ca 0 

Rah + ah} Raa~\~if aa I 


Case 3 
Case 2 
Case 1 


Z> f 0 1 rtl 

B?2 — “ 2^2 n 
R & 3 ^2 (^5 ~ r^o) 

^3 1 . 


5 , 




5/12 — 0 
Rh3 = 0 

C aX ^-a 

K 

Ca2 = 1 
Cc3 = 0 

c* 


^ 1 yi 

k 7 r 1 






(-*».) 


REAL EQUATIONS IMAGINARY EQUATIONS 


Rq aRafi 


1 aad a& 

R baR b& 


Ibaf b& 

R b&Rbh 

\ f bpf bh 

R e&R eh 

i f c&f ch 

RchRca 


1 chf ca 

RahRaa 


1 ahl aa 


RaaRafi 
1 baf b($ 

+ 

1 aal a& 
R ba R b{3 

Rb&Rbh 

+ 

Ib&Ibh 

i f c&f ch 

R c$R ch 

RchRca 

+ 

1 Chi ca 

1 ahl aa 

RahRaa 


= 0 Case 3 
-0 Case 2 


= 0 Case 1 


Note. — Terms with bars contain X; terms without bars do not contain X, 

The 9 quantities R aa , R a p, etc., refer to the real parts 
and the 9 quantities I aa , I a p, etc., to the imaginary 
parts of the coefficients of the 3 variables a, 0, and h 
in the 3 equations A , B, C on page 10. Denoting the 
coefficients of a, a, and a in the first equation by p f 
q, and r, 

B ” + ila.-l[-l> + iil v +r QJ ] 


C $ 3 — 0 
Ch\ — hi 

K 

Cm = 0 
Ch 3 = b 

♦The factor £ or ~ ^ not included in these constants. See the expressions for 
the /2’s and I’s on next page. 


which, separated in real and imaginary parts, gives 
the quantities R aa and I aa . Similarly, the remaining 
quantities R and / are obtained. They are all func- 
tions of k or C(k). The terms with bars R aat Rb&, 
and R ch are seen to be the only ones containing the 
[ unknown X. The quantities Q and X will be defined 
| shortly. The quantities R and I are given in the 
l following list: 
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*•■= -,1„, + U..Y + 1 2(i + «)[(* - a)0 - ^’] (1 , 

1 — -'O' + pA*l + (* + — 2 1 - 

/?,,* = A.,,+ * 2(^ + ~)<7 (3) 

lt ba =-li al - j. ’ - - a ) 6 ’ - J 7 ] ^ 

< /?»*=• If» + pIW: UA 

Rbh = ~ J, (l (6) 


ltd- -O. 27’,,, ^;/’l (8) 

/?r*-=-C\, + Si»A'-|2ff (it) 


/ “--i[ 2 ( a+ 5)IG- a > ,+ i°|-' 1 -] 

I a „ = — £ 2 (a + A' 


1 — A 
Ib “ k\ 


+*-] 

/ - 1 ! 
hf> ~k \ 

[SC 

1+®*] 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 



1 

~k 


G-“W"i 



1 

k 


J u F+->It, 0 (!) 

+c«] 


1..-1-2F 


(17) 

(IH) 

( 19 ) 


The solution as given by the three-row determinant 
shall be written explicitly in X. We are immediately 
able to put down for the general case a cubic equation 
in X with complex coeflicients and can easily segregate 
the three subcases. The quantity I) is as before the 
value of the determinant, but with the term containing 
X missing. The quantities J\f ua , M b t, and M ch are 
the minors of the elements in the diagonal squares 
aa, 6/3, and ch, respectively. They are expressed ex- 
plicitly in terms of R and / under the subcases treated 
in the following paragraphs. 


A an ~\~ “a A 



A,* 

A hfl -r%X 

Act 


A ah 
Abh 


= 0 


where A aa =R aa -\-iI ua etc. 


Comple: 

y cubic equation in 

X. 


UJWnX 

| A. 4 

L ria“P ■ ^h“-a A f,s) A 



r (12, i <-U/m 


(XXI) 

Case 3, 

(a, d) ■ 



SlailjX* r 

- ( A bt + -h A an) A - 

f -.(/,*=<! 

(XXII) 

Case 2, 

(0, h) 1 : 



Qsii.A’M 

- j- il fi Abt) X -) 

Af aa = i) 

(XXI11) 

Case 1, 

[ft , u): 



Q,,<K\'*+(<hA aa : -SJ.^A'H 

3/ ^ — 0 

(XXIV) 


£i " A k-\!rK ( 

w a /' a y 1 / br r o3 r \ 2 

itift r/ K \ t'k / 



' i<A “F.\/rV. = ( 

«£d\* 1 ( hr ' a '\ 

iOj-t’r) K \ t'k / 



v-JSei 

on, y 1 / br,u> T \ 2 
u r r T ) k \ r/r / 



and finally 



We are at liberty to introduce the reference param- 
eters (j) r and r ry and the convention adopted is: a> r is 
the last a) in cyclic order in each of the subcases 3, 2, 
and 1. 

Then 12„ = ( a — — ) and ll n +\ — 1, thus for 

\ w » \ \ r nj 1/ 

Case 3, 1 l a =(~ a '- a ) and 1 1 3 = 1 
(■use 2, and Si A =l 

Case 1, S2* = ^ w * ^ and il u — J 

To treat the general case of three degrees of freedom 
(equation (XXI)), it is observed that the real part 
of the equation is of third degree while the imaginary 
part furnishes an equation of second degree. The 
problem is to find values of A" satisfying both equa- 
tions. We shall adopt the following procedure: Plot 

graphically A” against | for both equations. The points 
of intersection are the solutions. We are only con- 
cerned with positive values of ^ and positive values of 

X. Observe that we do not have to solve for k, but 
may reverse the process by choosing a number of 
values of lc and solve for A". The plotting of X 

against for the second-degree equation is simple 

enough, whereas the task of course is somewhat more 
laborious for the third-degree equation. However, 
the general case is of less practical importance than 
are the three subcases. The equation simplifies con- 
siderably, becoming of second degree in X. 
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We shall now proceed to consider these three sub- 
cases. By virtue of the cyclic arrangement, we need 
only consider the first case (a, 0 ). The complex 
quadratic equations (XXII)- (XXIV) all resolve 
themselves into two independent statements, which 
we shall for convenience denote “Imaginary equa- j 
tion” and “Real equation”, the former being of first i 
and the latter of second degree in X. All constants ! 
are to be resolved into their real and imaginary parts, 
denoted by an upper index 11 or /, respectively. 

Let \f aa =M R aa -f iM r aa and let similar expressions j 
denote M and M ch 

Case 3 , (a, 0 ). 'Separating equation (XXII) we obtain. 

(1) Imaginary equation: 

(12 a 7*3 + ^ 3/41 a ) X + M T c h = 0 

v== _ ML* 

12 Q / / ) 3 -[-S 23 / a a 

( 2 ) Real ecjuation: 

il a Q$X 2 (il a Rbfi~\~ tt&Raa) A + M R eh =\) 

Eliminating X we get 

(Ujin+Ufilt'a) (Uaht+lkLJM 1 * 

+ M K a(iU b fi+n$Iaay = 0 

By the convention adopted we have in this case: 

Si " = (S)(S)’ and ^ =1 

Arranging the equation in powers of tt a we have: 

12/1 - A l' * (£*/*) + A/* rt V ] + «.[ {Mchf 
— M T c /i {Raal b$-\~ laaRbp) + 2 A/ 74 £*/„«/ 

+ [ ~ M r ( Jta«Iaa + - 0 

But we have 

~ - A/Ljj \R aa I b3 Ra&f 6a - ! - Rbfilt aa R bad "R aal bfi R b$I aa\ 

— M 1 ch {R a &I ba-V J a&Rb^) 

Finally, the equation for Case 3 (o, fl) becomes: 

12/ (M* ch I bf - A V Jlbtlbt) + «-[ ■ “ M' ch {Ra&ha + LfiRba) 
+ 2 Af^/ aa /^l + Af**/./-Af / rt ff^.=0 (XXV) 

where 

Af eft — RtJt* Ra&R'btt / aal bff 1 I afil ba 
A l ch — Raal b3 Rafi^ba~\ ' ' ,^R b 0 laftRba 

The remaining cases may be obtained by cyclic 
rearrangement : 


Case 2 , (ftA) r f 

V.Q{M R &a I c* — M l aa R C hl ch) 123 ! — Mi a (R (,*/ f/j-h I hhRcs) 

+2 MUbsIc^MLhy-MU^Jb^i) (A XVI) 

wheie A/ aa — R b ^Rtu RbhRcii I s/ca - ! - f j? 

A I a a ' R f> 3 ^ ch R bh I <*3 ~ I 1 ” / I bftR r3 

Case 1 , {hyOL) oc T = u) a 12 /I = ^ a, - N ) -- 2 12«=1 

\w a / r a 

I 2 jJ(*V/“^/ , aa 2 — H” ^ftl “ Aflg(R ea I a h ~T IcaRah) 

-f 2 A/« M /^/ a /-fd/^ 2 -d/L/^J cA = 0 (XXVII) 
where A f ^g — R C hR aa R ca Rah Ichl(ia~\" Icalah 

A I (,3 == R chi aa fC I ah ~F* / chRtia I caR ah 

Equations (XXV), (XXVI), and (XXVII) thus 
give the solutions of the cases: torsion-aileron, aileron- 
deflection, and deflection-torsion, respectively. The 
quantity 12 may immediately be plotted against 

^ for any value of the independent parameters. 

The coefficients in the equations are all given* in terms 
of R and I, which quantities have been defined above. 
Routine calculations and graphs giving 12 against 

^ are contained in Appendix I and Appendix II. 

Knowing related values of 12 and p X is immediately 

expressed as a function of 12 by means of the first- 
degree equation. The definition of A r and 12 for each 
subcase is given above. The cyclic arrangement of 
all quantities is very convenient as it permits identical 
treatment of the three subcases. 

It shall finally be repeated that the above solutions 
represent the border case of unstable equilibrium. 
The plot of X against 12 gives a boundary curve between 
the stable and the unstable regions in the A'S 2 plane. 

It is preferable, however, to plot the quantity p ^ 

instead of X, since this quantity is proportional to the 
square of the flutter speed. The stable area can easily 
be identified by inspection as it will contain the axis 

p V= 0 , if the combination is stable for zero velocity. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., May 2, 1934 • 
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APPENDIX I 

PROCEDURE IN SOLVING NUMERICAL EXAMPLES 


(1) Determine the li's and I's, nine of each for a 
major case of three degrees of freedom, or those per- 
taining to a particular subcase, 4 R’s and 4 1\ s\ Kefer 
to the following for the R’s and Fs involved in each 
case: 

The numerals 1 to 9 and II to 19 are used for con- 
venience. 

(Major case) Three 
degrees of freedom 

Raa La 11 

Ra? Lip 12 

Ran hn 13 

Rbti 1 ba 14 

R bp I bp 1 

R bti Lh 1 1> 

R(a La 17 

Re p Ice 1 8 

Ren 1 cn 19 

(Case .3) Torsional- 
aileron (a, 0) 

1 Raa Iaa H 

2 Rap he 12 

4 R ba Iba 14 

5 R b p I (,p lo 

(Case 2) Aileron - 
deflection (/?, h) 

5 R b p La 15 

G Rbh Lh Ih 

8 R c p he 18 

9 R ch I cb 19 

(Case 1) Deflection- 
torsion {h, a) 

7 Rea La 17 

9 R (h Lk 19 

1 R aa La 11 

3 Ran Ian 13 

It lias been found convenient to split the R’s in two 
parts R~R' ~rR" , the former being independent of 

the argument p The quantities I and R " are func- 


tions of the two independent parameters a ami r only. 6 
The formulas are given in the following list. 



U.i 

*'■<- f >•+ r„) p («+ [)(?'„«- 1 tJ) j 

(2) 

R "‘‘ h ~l 2 ( a + l) 6 ’ 

(?) 

R b ° k 7 r IV 2 ) ’ k j 

(4) 

/r M = ' k r- {■! ( T " a ~i TwF ) ~ fc 1 Tr ~ Ti 7 ’ lo) } 

(5) 

, r 1 T a 

K nn — T --- 

fi 7 r 

(6) 


(7) 

R''^ ll(r u a- 2 r u ?) 

(»; 


(») 

/„„ = — 2 (a + g) i( 2 - «) f + J- 0 } t 2 - " 

(ii) 

I „ a = — “ |( « + k 

(12) 

+Q «) r ] 


I ah = 2 (a + 2 ) F 

(13) 

Whore p — -- J (1 - c*) al * 

(14) 

h>=L L( TuF+ l r 'L 

(15) 

/«-— r 

7T 

(16) 


(17) 

e-T 

1 

O 

<m :-4e 
+ 

ii 

HQ 

<- 

(18) 

I eh ~ 2 F 

(19) 


— s <yhe quantities / given in the appendix and used in the following calculations 
are seen to differ from the Is given in the body of the paper by the factor j, • It 
may be noticed that this factor drops out in the first-degree equations. 


1 

2 

3 

4 

5 
G 

7 

8 
9 


302 



GENERAL THEORY OF AERODYNAMIC INSTABILITY AND THE MECHANISM OF FLUTTER 


Choosing certain values of a and c and employing 
the values of the T’s given by the formulas of the report 
(p. 5) or in table I and also using the values of F and 
G (formulas (XII) and (XIII)) or table II, we evaluate 

the quantities I and R" for a certain number of j- 

v allies. The results of this evaluation are given in 
tables III and IV, which have been worked out for 
a = 0,— 0.2, and — 0.4, and for c=0.5 and c — 0. The 

range of is from 0 to 40. These tables save the work 

of calculating the I* ft and R"'s for almost all cases of 
practical importance. Interpolation may be used for 
intermediate values. This leaves the quantities R' to 

be determined. These, being independent of are as 

a result easy to obtain. Their values, using the same 
system of numbers for identification, and referring to 
the definition of the original independent variables on 
pages 9 and 10, are given as follows: 




(i) 

ir ma — r £-{e-a)« + 

K K 

J j + (c- «)" 

(2) 

It'an a 

K 


(3) 

R' ba = same as R' flt s 


(4) 

v + > 


(5) 

H\h= + 

K 7 T 


(«) 

R' ra — same as R' fth 


(7) 

/£+=same as R\ h 


(8) 

w, »=---l 

K 


(9) 

‘ of the symmetrical 

arrangement in 

the 


determinant, the 9 quantities, are seen to reduce to 
0 quantities to be calculated. It is very fortunate, 
indeed, that all the remaining variables, segregate them- 
selves in the f> values of R' which are independent of y 


while the more complicated / and R" are functions 
solely of c and a. In order to solve any problem it is 
therefore only necessary to refer to -tables III and IV 
and then to calculate the 6 values of R f . 

The quantities (1) to (9) and (11) to (19) thus 

having been determined, the plot of 12 against p which 

constitutes our method of solution, is obtained by 
solving the equation aU 2 +bil+c~ 0. The constants 
a, 6, and c are obtained automatically by computation 
according to the following scheme: 


Case 3 


Find products 1.5 

2.4 

11,15 

12.14 

Then M R ch = 1 . 5 — 2 .4 

— p(H.15- 

-12.14) 


Find products 1.15 

2.14 

11.5 

12.4 


Then M 7 cA =l. 15-2. 14 + 11.5-12.4 
and a=M R ch(lo) 2 — M 7 cA (5.15) 

6 = -M 7 , A (2.14 + 12.4)+M+ A (11.15) 
c=M R ch ( 1 1 ) 2 — A/ 7 cA ( 1.11) Find 12* 

Solution: 

A ill cfl 

Similarly 

Case 2 

5.9 0.8 15.19 16.18 

M*,«=5.9— 6.8 — p(15. 19— 10.18) 

5.19 6.18 15.9 16.8 

A/ 7 ,,— 5.19— 6.18+15.9— 16.8 
a = M“ aa ( 1 9) 2 — M' a „ (9.19) 
ft = — il/ 7 „„(6.18— 16.8) — 2A/ 7 '„ 0 (15.19) 
c = Af “ a «(l 5) 2 —M‘ aa (5.15) Fiiyi U 3 

1 .0,(19) + 15 

X ~ M’a . 

and 

Case 1 

9.1 7.3 19.11 17.13 

M>=9.1-7.3-p(19.1 1-17.13) 

9.11 7.13 19.1 17.3 

A/ 7 # 0=9.1 1 — 7.13 + 19.1 — 17.3 
a=Af R „0(ll) 2 — A/'^l.ll) 
b - - AF „0(7. 13 + 17.3) + A/“„(19. 1 1) 
c = A/ R 40 ( 1 9 ) 2 — A/ 7 j0 (9.19) Find Sl„ 

1 Q»(ll) + 10 

:v m‘ 00 

<X is defined as f for case 3 ; 

lie is defined as for case 2 > and 

12 ft is defined as for case 1 . 

The quantity ^ is «(^r) b 3 T definition. 

Since both 12 and are calculated for each value of 

A, we may plot ^ directly as a function of 12. This 

quantity, which is proportional to the square of the 

flutter speed, represents the solution. 

We shall sometimes use the square root of the above 

quantity, via, J will denote thi. 
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quantity by F, which we shall term the “flutter factor ” 
The flutter velocity is consequently obtained as 


v = F 


bu T r r 

v« 


ba) T r, 


Since F is nondirnensional, the quantity must 

V«f 

obviously be a velocity. It is useful to establish the 
significance of this velocity, with reference to which 
the flutter speed, so to speak, is measured. Observing 

that k = and that the stillness in case 1 is given by 

/ -Tjrf--* this reference velocity mav be written: 

V Mb~rJ 

_bu lt r a i jc 
l R~ — ,1 — or 

\K 0 \ 7 rp 


it pty*6 2 = C,, 

The velocity is thus the velocity at which the total 

force on the airfoil irpv^Zb attacking with an arm ^ 

equals the torsional stiffness C a of the wing. This 
statement means, in case 1, that the reference velocity 
used is equal to the “divergence” velocity obtained 
with the torsional axis in the middle of the chord. This 
velocity is considerably smaller than the usual diver- 
gence velocity, which may be expressed as 


1 

Vi>**v R i 

2 + a 

where a ranges from 0 to - We may thus express 
the flutter velocity as 

V F = 


In case 3 the reference velocity has a similar signifi- 
cance, that is, it is the velocity at which the entire lift of 

the airfoil attacking with a leverage ~ b equals numeri- 
cally the torsional stiffness C$ of the aileron or movable 
tail surface. 

In case 2, no suitable or useful significance of the 
reference velocity is available. 


TABLE I.— VALUES OF T 


* c - 1 

c='A 

C=0 

c= — l A 

c=-l 

T: 

Tj.- 

1 0 

0 

-0. 12.59 
-0.2103 

-0. 6667 
-1.5707 

-1.6967 
-4. 8356 

-3. 1416 
-9. 8697 

Ti ... 

0 

05313 

8084 

-3.8375 

-11. 1034 

74.. 

. . . . . 0 

-.6142 

-1. 5 70S 

-2. 5274 

-3. 1416 

T S 

: 0 

— . 9398 

-3. 4674 

-6. 9.503 

-9. 8697 

n 

J 0 

-0 2103 

— 1. 5707 

-4. 8356 

! -9. H697 

T; 

! o 

.0132 

1964 

-1. 1913 

-3. 5343 

1 n . . 

; (J 

.0903 

-. 3333 

- 1. 480.5 

-3. 1416 

I 7’in ... 

0 

3.9132 

2.5708 

2.9604 

3. 1416 

U. 

(J 

1. 2990 

3. 570S 

6. 3538 

9. 4248 

! Tu . 

0 1 

. 07066 

4292 

1. 2990 

3. 1416 


TABLE II. — TABLE OF THE BESSEL FUNCTIONS J tu J { 
>0. Y t AND THE FUNCTIONS F AND G 




k 

1 

k 

./« 

Ji 

Ya 


F 

—a 

to 

0 





0. 5000 

o 

30 


- 0. 2459 

0. 0435 

0. 0.557 

0. 2490 

. .5006 

0.0126 

6 

u 

. 1501 

-.2767 

-.2882 

1750 

. 5018 

.0207 

4 

u 

- . 3972 

-.0660 

-.0170 

. 3979 

. 5037 

.0305 

2 

!»'/ 

2239 

. 5767 

.5104 

1071 

.5129 

.0577 

I 

i 

. 7652 

. 4401 

.0882 

7813 

.5395 

. 1003 

.8 

m 

. Mf>3 

.3088 

— . 0Wi8 

9780 

5541 

. 1165 

. 6 


.9120 

.2867 

-. 308.5 

-1. 2604 

. 5788 

. 1378 

. 5 ' 

2 

. 9:485 

.242:4 

-. 4444 

-1. 4714 

. 0030 

. 151 

. 4 

2j>$ 

.9604 

. I960 

- . 0060 

- 1. 7808 

. 6245 

. 166 

• » 

m 

. 9776 

. 1483 

- . 8072 

-2. 2929 

. 6650 

. 180 

. 2 

5 

. 9900 

0995 

— 1. OSH) 

-3. 3235 

. 7276 

. 1886 

■ 1 

30 

. 9975 

. 0499 

-1-5342 

- 7.0317 1 

. 8457 

. 1626 

. 05 

20 





.911 

1 . 132 

.025 

40 





. 965 

| .090 

0 

oo 





1.000 

: 0 
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TABLE TIL— VALUES OF R 


k 

0 

M0 


M 

W 

1 

I 

1 M 

m 


m 


5 

10 

20 

40 

j 

c 

a 







1 




0 

0 

-0.00564' 

-0.01566 

-0. 03529 

-0. 14265 

(-0.58965 

-0. 93656 

-1.72330 

-2. 56300 

-4. 11000 

— 7. 68720 

! 

-18 66150 

-85. 38300 

-365. 72000 

-1, 528. 2000 

R" aa 

(') 

-.2 

0 

- 00353 

-.00981 

-.02208 

-.08905 

: -.36586 

-.58061 

- 1.061.58 

-1.5740'- 

-2.51580! -4.68430 

-11.31010 

- 51. 4241*0 

-219. 74900 

-917.3520 


. 4 

0 

001-23 

-.00341 

00767 

j - .03084 

, -. 12595 

19936 

-.36305 

53676 

-.85520 -1.58540 

-3. 80774 

j -17.20670 

-73. 35520 

-305. 9280 



o| 

0 

-.001631 

-.00452 

01020' 

-.04175 

' 18016 

-.29384 

-.56223 

-.87212 

-1. 43983 

-2. 84988 

-7. 46300 

-38.29650 

-172.36360 

-741. 7972 


0 

-.2 

0 

. <xxxw 

. 00083 

.00184 

. 00679 

. 01922 

. 02266 

.01629 

— .01 400 

-.06.803' -.29517 

- 1. 29480 

| -10.24590 

-52. 49020 

-241.3664 

R'hd 


-.4| 

0 

. 00-222 

.00617 

j . 01388 

05531 

.21861 

.33914 

. 59499 

.84414 

1. 30365 

2. 25914 

4. 87340 

17.80470 

67. 38320 

j 259. 0648 


0 

0 

. 00083 

. 00229 

1 . 00510 

. 01932 

. 06419! 

. 08876 

. 12176 

. 12260 

12205 

- . 02900 

-.93535 

-10. 48970 

-59 16180 

-268. 7236 


0.5 

-.2 

1 0 

.00214 

.00595! .01336 

. 05278 

. 20325! 

.31065: .53062 

. 73222 

1. 10233 

1.81135 

3. 55230 

10. 14740 

. 31.49620 

101. 6340 



4 

0 

.00347 

. 00965 

.02170 

. 08656 

, 3436 1| 

| . 53463 

! . 94336 

1.34762 

2. 09190 

3.66913 

8. 08235 

30. 97980 

120. 89760 

476. 2592 



0 

0 

-.00125 

-.00345 

: -.00763 

-.02890 

-. 10030 

j -.145601 -.22470 

- 30200 

-.41500 

-.60000 

-.94300 

: -1.62600! — 2. 64000 

-3. 6000 


(') 

-.2 

0 

-.00075 

-.00‘207 

- . 00426 

| -.01734 

. 1)6018 

- . OH736 

-. 13482 

18120 

-.24900 

-.36000 

-.56680 -.97560 

-1.58400 

-2. 1600 



- 4 I 

0 

-.00201 

-.00334 

-.00.502 

-.01003 

-.02006 

-. 02508 

! -.03236 

-.04012 

-.05015 

-.06683 

-. 10030 

20060 

-. 40120 

-.8024 



o! o 

.00077 

.00214 

.00482 

. 01949 

. 08055 

. 12821 

.23541 

.35010 

56143 

1. 05008 

2. 54920 

11.66330 

49. 95700 

1 208. 7520 


0 

- 2 

O' 

. 00080 

. 00223 

. 00503 

. 02027 

. 08329 

. 13219 

.24169 

. 35836 

. 57276 

1 . 06650 

2, 57490 

11. 70770 

50. 03000 

208.8500 

R"t a 


4 

| 0 

.00084 

. 00233, 

. 00523 

. 02106 

. 08603 

. 13616 

.24796 

. 36661 

.58410 

1.08286 

2. 60069 

11.75220 

50. 10160 

208. 9490 


(J 

0 

. 00013 

. 0003 5 1 

.00079 

.00321 

. 01327 

.02112 

. 03878 

. 05767 

. 09248 

. 17296 

. 41988 

1. 92110' 

8.22870 

34. 3850 


0.5 

-.2 

0 

.00013 

. 00037 

.00083 

. 00334 

. 01372 

.02177 

. 03981 

.05903 

. 09434 

. 1 7 566 ' 

. 42413 

1.92840 

8. 24060! 

34. 4007 


1 

4 

0 

. 00014 

.00038 

. 00086 

. 00347 

.01417 

! . 02243 

. U40S4 

.06039 

.09621 

. 17836 

. 42837 

1. 93575 

8. 252461 

34. 4369 

R” b a 

i o 


0 

. 0012 4 

' . 00343 

.00772 

.03101 

. 12642 

. 19830 

. 35807 

. 52400 

. 82930 

1. 5163 

3. 54970 

15 351201 

64. 02240 

263. 2340 


<*) 

1} 

.00031 

. 00087 

.00196 

. 00785 

.03170 

. 04980 

. 08935 

. 13000 

20440 

. 3694(1 

. 84970 

3. 55050 

14.56740 

59. 3188 

r>" 

0 


0 

.00017 

1 .00047 

.00104 

.00394 

.01370 

. 01989 

. 03177 

.04125 

. 05669 

. 08196 

. 12881 

. 22211 

. 36062 

. 4918 

bk 

, 5 

(>) 

0 

.00003 

.00008, 

. 00016 

f . 00065 

. 00220 

. 00328 

.00506 

.00680 

. 00934 

01350 

.02122 

03659 

. 05940 

.0810 



0 

0 

.01)28| .03)321 

. 07058 

| . 28530 

i 1. 17930 

1.87710 

3. 4467C 

5. 12600 

8. 22000 

15, 37450 

37. 32300 

170. 7 6600 

731. 44000 

3, 066. 4000 

R"ta 

(■) 

-.2 

c 

.01176' .03270: 

. 07362 

. 29684 

! 1. 21954 

1. 93540 

3. 53860 

5. 24680 

8. 38600 
H. 55200 

15.61440 

37. 70020 

171. 41640 

732. 49600 

3,057.8400 



4 

0 

.01226; 

. 03408 

. 07668 

. 30838 

1. 25950 

! 1.99360 

3. 63050 

5. 36760. 

15.85440 

38. 07740 

172. 06680 

733. 55200 

3, 059. 2800 

R"tfi j 

0 

.5 

(*) 

0 

0 

. 00003 
. 00660 

. 02673 
.01840 

. 06018 
.04150 

. 24266 
. 16810 

~~ 1.00561 
. 69850 

1. 58£)6 
1. 11453 

2. 89371 
2. 0532(1 

4. 29100 
! 3. 06224: 

6. 85898 
4. 92530 

1549965 
9. 24438 

30. 84330 
22. 51400 

140. 26370 
103. 67300 

599. 41300 
444. 86400 

2, 502! 3470 
1,881.990(1 

.fjiJ 


( 3 ) 

0 

. 00250 

. 00690 

.0142(1 

. 05780 

. 20060 

. 29120 

. 44940 

| . 604001 

. 83000 

1. 2(XXX1 

1.88600 

3. 25200; 

5. 28ypo 

~~ 7~200U 


1 Independent of c. 


1 Independent of a. 


TABLE IV— VALUES OF I 


1 

k 

0 

Mo 


H 



1 M 

1 ^ 

2 

2 M 

3M 

5 

10 

20 

40 


c 

a 






0 

0. 25000 

0. 2,5096 

0. 25255 

0. 25578 

0. 27240 

0. 33055 

0.36S55 

0. 44030 

0. 50050 

0. 60275 

0. 76750 

1.07920 

1. 70320 

2. 68450 

3.61750 

Ia a 

(>) 

-.2 

. 49000 

. 49050 

. 49131 

. 49302 

. 50189 

. 53359 

. 55464 

. 59472 

. 62794 

. 68671 

.78070 

. 96021 

1.32040 

1.90140 

2.45470 



-. 4 

.81000 

.81014 

. 81037 

.81086 

.81145 

.82395 

. 82938 

.84176 

.85186 

. 87059 

.90030 

. 95763 

1. 07300 

1. 26400 

1.44630 



0 

. 17805 

. 17874 

. 17985 

. 18219 

. 19433 

. 23768 

. 26645 

.32132 

.36664 

. 44690 

. 57526 

. 82035 

1.31213 

2. 10476 

2. 85963 


0 

-.2 

39170 

.39212 : 

. 39278 j 

.39418 

.40147 ! 

. 42748 

. 44474 

.47761 

. 50485 

.55300 

. 63002 

. 77708 

1.07215 

1. 54773 

2. 00065 

lad 


- 4 

.60531 

.60545 , 

.60567 

.60614 

. 60857 

.61724 

. 62299 

.63395 

. 64303 

. 65908 

. 68475 

.73377 

.82313 

. 99065 

1. 14163 


0 

! . 13252 

. 13317 

. 13425 

. 13640 

. 14742 

.18544 

. 20914 

.25611 | 

. 29514 

. 35951 

. 46379 

.65973 

! 1.05124 

1. 65524 

2. 22869 


0.5 

-.2 

I .21297 

. 21336 

.21401 

. 21530 

.22191 

. 24472 

i . 2.5894 

.28712 

.31054 

. 34916 

.41173 

. 52929 

. 76420 

1. 12651 

1.47067 



-.4 

. 29342 

.29354 

. 29376 1 

. 29419 

. 29640 

.30400 

.30891 

.31813 

. 32594 

. 33881 

. 35966 

. 39884 

.47714 

.50792 

. 71260 



0 

1 — . 50000 

- . 50060 j 

!-. 50180 j 

-.50370 

-.51290 

-.539.50 

-.5.5410 

-.57880 

i -.60300 

-. 62450 

-.66500 

72760 

-. 84570 

-.94100 

-.96500 

/.* 

(■) 

-. 2 ; 

- . 30000 

- . 30036 j 

30108 1 

-. 30222 

-.30774 

-.32370 

-.33246 

-.34728 

; -.36180 

37470 

-.39900 

43656 

50752 

-.56460 

- . 57900 



-.4 

10000 

-. 10012 1 

! — . 10036 

-. 10074 

10258 

-. 10790 

1 1082 

11576 

| 12060 

12490 

-. 13300 

14552 

-. 16914 

-. 18220 

19300 



0 

i .39023 1 

.39010 

. 38988 ! 

. 38944 

.38717 

. 37923 

. 37404 

. 36424 

. 35601 

.34204 

.31954 

> . 27696 

. 19172 

. 05766 

-.06980 


0 

-.2 

. 40389 

. 40378 ! 

. 40359 

. 40320 

.40119 

. 39397 

, . 38918 

.38005 

. 37249 

. 3591 1 

.33771 

. 29683 

.21469 

. 08255 

-. 04344 

/bo 


-.4 

.41755 

.41746 

. 41730 

41697 

.41520 

.40871 

j . 40432 

. 39586 

.38896 

.37fd7 

. 35598 

.31671 

. 23793 

. 10744 

-.01707 


0 

. 07438 

. 07435 

. 07433 

.07424 

.07387 

.07256 

.07171 

.07009 

.06874 

.06644 

.06273 

. 05572 

.04168 

.01960 

-.00327 


0.5 

-.2 

. 07663 

. 07661 

. 0765ft 

. 07651 

.07618 

.07499 

. 07420 

. 07270 

.07145 

.06925 

. 06572 

. 05899 

.04548 

.02370 

.00295 



- 4 

07887 ' 

; . U7885 

. 07882 

j .07867 | 

.07848 

.07741 

. 07668 

.07529 

.07416 

. 07205 

.06871 

. 06226 

.04928 

. 02779 

.00728 

Ufi 

0 

<*) 

.32297 i 

. 32288 

. 32273 

.32241 

. . 32075 

.31483 

.31090 

.30342 

. 29721 

. 28625 

. 26872 

. 23524 

. 16806 

. 05979 

-.04333 

.5 

.04270 

. 04270 

. 04270 

. 04270 

,04240 

.041,50 

. 04095 

. 03930 

. 03904 

.03760 

.03386 

.03080 

.02200 

.00845 

-.00470 

/»A 

0 

(’) 

. 06830 

.06840 

. 06850 

.06880 ; 

. 07010 

. 07370 

. 07570 

.07910 

. 08240 

.08530 

.09080 

.09940 

.11550 

. 12440 

. 13180 

. 5 

| .01125 I 

.01126 

.01129 

.01133 1 

. 01 154 

.01214 

.01247 

.01302 

.01357 

. 01405 

.01496 

.01637 

01903 

.02117 

.02171 



0 

1. 50000 

1. 49808 

1.49490 

1. 48844 

1. 45520 

1 33890 

! 1. 26290 

1. 1194(1 

.99900 

. 74950 

. 46500 

-. 15840 

-1 40630 

-3. 36900 

-5. 23500 

ha 

(’) 

-.2 

1. 70000 

1 . 69832 

1.69562 

1.68992 

1.66036 

1.55470 

: 1. 48454 

1.35092 

1 . 24020 

1.04430 

. 73100 

. 132G4 

-1.06802 

-3. 00460 

-4.84900 


-. 4 

1. 90000 

1. 89856 

1.89634 

1.89140 

1.88552 

1 77050 

! 1.70618 

1.58244 

1.48410 

1. 31410 

.99700 

.42370 

- 72974 

-2.64020 

-4.46300 

he 

0 

(*) 

I. 06830 

1. 06090 

1.06470 

1.06000 

1.03580 

.94900 

.89150 

.78190 

. 691 10 

. 52840 

. 27380 

-.21640 

- 1 . 20010 

-2. 78550 

-4.29530 

.5 

. 40220 

.40100 

. 39880 

.39450 

. 37240 

.29640 

. 24640 

. 15510 

.07610 

-.05180 

26030 

-.65220 

-1.43520 

-2.64380 

-3. 79010 

Lk 

0) 

( 2 ) 

1.00000 

1.00120 

1.00360 

1.00740 

1.02580 

1.07900 

1. 10820 

| 1. 15760 

1 . 20600 

1. 24900 

1.33000 

1.45520 

1.69140 

1.82200 

1.93000 


1 Independent of c. 


1 Independent of a. 



APPENDIX II 


NUMERICAL CALCULATIONS 


A number of routine examples have been worked out 
to illustrate typical results. A “standard” case has 
been chosen, represented by the following constants: 


k = 0. 1 , c=0.5, a= — 0.4, x a = 0.2, 


1 


»•«*= 0.25, = ^£> 2 = infT 


100 


u> a , 03/3, variable. 

We will show the results of a numerical computation 
of the three possible subcases in succession. 



Case 3, Torsion-aileron (a, 0 ): Figure 5 shows the <J a 

against ^ relation and figure t> the final curve 

F=k(- 17 vYagainst Wa -V=4o(°^Y 

\c*) 3 r s by Kuprt}/ \a)<}/ 



Figure fi.— Case 3. Torsion-aileron (or. (J): Standard case. Showing flutter factor 
F against U.>. 


Case 2 , Aileron- flexure ($, A): Figure 7 shows the 
against relation 6 and figure <S the final curve 



against 

\ / 



4 It is realized that considerable care must he exercised to get these curves reason - 
bly accurate. 


The heavy line shows the standard case, while the 
remaining curves show the effect of a change in the 

value of Xr to 77 and 7777 
40 100 

Case 1 , Flexure-torsion (A, a): Figure 9 shows again 



Figure 7. -Vase 2. Aileron-deflection (fl , h) . (a) Standard case, (b), (c). (d) indicate 
dependency on it j. Case (d), x& = —0.004, reduces to a point. 


the S2„ against 


1 

k 


relation and figure 10 the final result 



against = -*(£) 


Case 1 , which is of importance in the propeller theory, 
lias been treated in more detail. The quantity F shown 


in the figures is \k 


Figure 11 shows the dependency on^=- 1 ; 

figure 12 shows the dependency on the location of the 
axis a; figure 13 shows the dependency on the radius of 
gyration r a =r; and figure 14 shows the dependency 
on the location of the center of gravity x } for three 
different combinations of constants. 

EXPERIMENTAL RESULTS 

Detailed discussion of the experimental work will not 
be given in this paper, hut shall be reserved for u later 
report. The experiments given in the following are 
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restricted to wings of a large aspect ratio, arranged with I able springs restrain the wing to its equilibrium 


two or three degrees of freedom in accordance with the 



fL a 

Figure 8. —Case 2, Aileron-deflection (&, A): Final curves giving flutter factor h 
against lift corresponding to cases shown in figure 7. 

theoretical cases. The wing is free to move parallel to 
itself in a vertical direction (h); is equipped with an 



0 2 4 6 8/0 

i/k 

Figure ».-Oiise 1, Flexure-torsion (A. «): Standard case. Showing ilk against y 

axis in roller hearings at (a) (fig. 2) for torsion, and 
with an aileron hinged at (c). Variable or excliange- 


position. 



Figure 10.— Case 1, Flexure-torsion (A, «): Standard case. Showing flutter factor 
F against Si*. 

We shall present results obtained on two wings, both 
of symmetrical cross section 12 percent thick, and with 
chord 26—12.7 cm, tested at 0°. 



Figure II. -Case l, Flexure-torsion (A, «): Showing dependency of F on^* The 
upper curve is experimental. Airfoil with r— ^ a = -0-4; r- 0.2; 4« =. 01, — variable. 

Wing A, aluminum, with the following constants: 

1 - 0.4, jv-=0.31, 0.173, and 0.038, 

410 

respectively ; 

rJk-0.33 and u* = 7X2ir 
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Wing B, wood, with flap, and the constants: 

k = c= 0.5, a= — 0.4, z„=0.192, r a 2 = 0.178, 

*0=0.019, r a 2 = 0.0079, and kept constant 

= 17.6X2x 

The results for wing A, case 1, are given in figure 15; 
and those for wing B, cases 2 and 3, are given in figures 
16 and 17, respectively. The abscissas are the fre- 
quency ratios and the ordinates are the velocities in 
cm/sec. Compared with the theoretical results calcu- 
lated for the three test cases, there is an almost perfect 
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Figure 12.— Case 1. Flexure-torsion ih, «): Showing dependency of F on location 

of axis of rotation a. Airfoil with r= x = 0 2 * - — h; a Variable 
2 4 «2 6 

agreement in case 1 (fig. 15). Not only is the minimum 
velocity found near the same frequency ratio, but the 
experimental and theoretical values are, furthermore, 
very nearly alike. Very important is also the fact that 
the peculiar shape of the response curve in case 2, pre- 
dicted by the theory, repeats itself experimentally. 
The theory predicts a range of instabilities extending 
from a small value of the velocity to a definite upper 
limit. It was very gratifying to observe that the upper- 
branch of the curve not only existed but that it was 
remarkably definite. A small increase in speed near 
this upper limit would suffice to change the condition 
from violent flutter to complete rest, no range of transi- 
tion being observed. The experimental cases 2 and 3 
are compared with theoretical results given by the 
dotted lines in both figures (figs. 16 and 17). 


The conclusion from the experiments is briefly that 
the general shapes of the predicted response curves re- 



Figure 13 —Case 1, Flexure-torsion (A, a): Showing dependency of Ton the radius 
of gyration r n = r. 

A, airfoil with a= — 0 . 4 ; j=o.2 : T variable. 



O .2 .4 .6 


Figure 14 —Case l. Flexure-torsion (A, a): Showing dependency of F on x a , the 
location of the center of gravity. 

A. airfoil with r = -!-V=-0.t; « -- d w . ' — I— -I • z variable. 

-s 4UU u>3 ■ 0 

B. airfoil with r-h; a = - 0.4: *= : — — 1 : * variable. 

i 4 u i 6 

C. airfoil with r-1: <, = -0 4; ,= A; ^=1; z variable, 

peat themselves satisfactorily. Next, that the influ- 
ence of the internal friction 7 obviously is quite appreci- 

* This matter is the subject of a paper now iu preparation. 
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able in case 3. This could have been expected since 
the predicted velocities and thus also the air forces on 
the aileron are very low, and no steps were taken to 
eliminate the friction in the hinge. The outline of the 
stable region is rather vague, and the wing is subject 



velocity r against frequency ratio Deflection-torsion. 

to temporary vibrations at much lower speeds than 
that at which the violent flutter starts. The above 
experiments are seen to refer to cases of exaggerated 
unbalance, and therefore of low flutter speeds. It is 
evident that the internal friction is less important at 
larger velocities. The friction does in all cases increase 
the speed at which flutter starts. 



Figure 16.— Case 2. Wing B. Theoretical and expert mental curves giving flutter 
velocity v against frequency ratio Alleron-deflentlon fji> A). 



Figure 17.— Case 3. Theoretical curve giving flutter velocity against the fre- 
quency ratio The experimental unstable area Is indefinite due to the im- 
portance of internal friction at very small velocities. Torsion-aileron (a, 0) 
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